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ABSTRACT 

We observed the double neutron star binary (DNSB) containing PSR 
J1537+1155 (also known as B1534+12) with the Chandra X-ray Observatory. 
This is one of the two DNSBs detected in X-rays and the only one where a hint 
of variability with orbital phase was found (in the previous Chandra observa- 
tion). Our follow-up observation supports the earlier result: the distribution of 
photon arrival times with orbital phase again shows a deficit around apastron. 
The significance of the deficit in the combined dataset exceeds 99%. Such an 
orbital light-curve suggests that the X-ray emission is seen only when neutron 
star B passes through the equatorial pulsar wind of neutron star A. We describe 
statistical tests we used to determine the significance of the deficit, and conclu- 
sions that can be drawn from its existence, such as interaction of the pulsar wind 
with the neutron star companion. We also provide better constrained spectral 
model parameters obtained from the joint spectral fits to the data from both 
observations. A power-law successfully fits the data, with best-fit photon index 
T = 3.1 ±0.4 and unabsorbed flux / = (3.2 ±0.8) x 10~ 15 ergs^cm- 2 (0.3-8 keV 
range) . 

Subject headings: binaries: close — pulsars: individual (PSR J1537±1155 = PSR 
B1534±12) — stars: neutron — stars: winds, outflows 



1. Introduction 

Double neutron star binaries (DNSBs) are rare systems comprised of two neutron stars 
(NSs) in a tight orbit. In addition to the gravitational interaction, in which general rela- 
tivistic effects often play a significant role, the two NSs may also interact in other ways (e.g., 
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via electromagnetic torque or pulsar wind). The manifestations of such interaction(s) can be 
found via precise radio timing, provided that one of the NSs is seen as a radio puls ar. In the 
exceptional case of DNSB J0737— 3039 (hereafter J0737; see lKramer fc Stairsll2008l for recent 
review), both the NSs are puls ars, which has enabl ed a suite of high-precision measurements 
of post-Keplerian parameters ( iKramer et al.ll2006l ). determination of the component masses 
( iLyne et al.ll2004j ). and anal ysis of the inte raction between the wind of pulsar A and the 
magnetosphere of pulsar B (ILyutikovl I2004T). The evolutiona ry processes in which DNSBs 
may form have been discussed by, e.g., Ivan den Heuvell ( 120051 ). 



J0737 had been t he only DNSB detected in X-rays (McLaughlin et all 120041 ) . until 
Kargaltsev et al.l (120061 hereafter KPG06) discovered X-ray emission from DNSB J1537+1155 
(hereafter J1537). KPG06 observed J1537 with the Advanced CCD Imaging Spectrometer 
(ACIS) on board Chandra. In the 36 ks exposure, they detected 16 photons, whose energies 
suggested a steep spectrum (photon index r « 3). Intriguingly, KPG06 found that the 
arrival times did not appear to be evenly distributed in orbital phase, showing a deficit of 
photons near apastron. However, the statistical significance of the deficit in the first ob- 
servation was marginal (in fact, as we show below, it was overestimated by KPG06). Here 
we present analysis of the second Chandra ACIS observation (with a re-analysis of the first 
one), which provides additional support for the dependence of its emission on orbital phase. 



In Section 2, we describe the data analysis, with details about the timing 
spectrum (§2.2). We present conclusions and a brief summary in Section 3. 



52.1) and the 



2. Observation and Analysis 



In both observations of J1537, the ACIS was operated in the Very Faint (VF) m ode. 
For each of the datasets, we applied VF-mode background filtering (IGarmire et al.ll2003l ). re- 
moved positional randomi zation applied by the pipeline, and app lied the sub-pixel resolution 
tool to split-pixel events (ITsunemi et al.ll2001t iMori et al.ll200lf ). 



The first observation (ObsID 5591) was carried out on 2005 April 10, with a total 
exposure of 36080 s (after filtering for good time intervals). The source was imaged on-axis on 
the most sensitive S3 chip, operated in Timed Exposure mode with a frame time of 3.24104 s 
(full-chip readout, giving 457s dead-time and 35623s net live time). We reprocessed this 
data set and generated a new Level 2 event fild^, in order to ensure consistency with a more 
recent pipeline version (DS 8.2) used to process the second data set. As a result of slightly 
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different processing, the extraction aperture of r = 1'.'23 (used in KPG06) now includes 17 
photons (one extra photon compared to KPG06) in the energy range of 0.3-10 keV. For 
spectral analysis, photons were extracted using the CIAO task psextract (CIAO version 
4.2, CALDB version 4.2.0). Counting the number of background counts in an annulus 
(4'.'92 < r < 44'.'28) centered on J1537, we expect to find (on average) 0.37 background 
photons within the source extraction aperture. 

The second observation (ObsID 20224) was carried out on 2009 September 27 and had 
a total exposure time of 38181 s. There were no bad time intervals. Again, the ACIS S3 chip 
was used, this time in 128 row sub-array mode (i.e., 1/8 of the chip), providing 441.04ms 
time resolution. The dead-time for such short frame time is 9.3%, so the total live time 
after the dead-time correction is only 34628 s. The less efficient high time resolution mode 
was used in order to make it possible to detect pulsations in the data, should slow radio 
pulsations from neutron star B be detected in the radio. As above, we extracted photons 
from an r = 1'.'23 aperture centered on the source. The total number of photons is 16, in the 
0.3-10 keV energy range. Of these, on average, 0.69 photons are expected to come from the 
background. The image from the second observation is shown in Figure [TJ There appears 
to be no detector artifacts in the vicinity of the source, and the background beyond the 
extraction circle is quite uniform. 



2.1. Timing 

For the timing analysis, we used the photon arrival times from the Level 2 event files, 
after correcting them to the solar system barycenter with the axbary tool from CIAO 4.2. 
For each photon, we calculated the orbital phase using the radio ephemeris (courtesy of 
Ingrid Stairs). The phase of each detected photon is plotted in Figure [2j There appears to 
be a portion of orbital phase in which few/no photons are detected from the source. We 
have checked that the fluctuations in the background rate were consistent with Poissonian 
statistics, and no decrease in the background rate was seen in the phase intervals with no 
counts detected in the source aperture. 

As an alternative representation of the same data, Figure [3] shows the phase-folded 
light curve, with three bins, where we choose one of the bins to be centered at the phase of 
apastron, 0^ = 0.5. We also show confidence intervals for each bin. 

We will now assess the probability that the observed deficit of photons is caused by 
random fluctuations, assuming that the number of counts within a certain fixed time interval 
obeys Poissonian statistics. 
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Fig. 1. — 40"x40" image of the vicinity of J1537 in the 0.3-10 keV range, obtained in the 
second observation. The location of the r = l'/23 extraction circle is shown in green. 
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Fig. 2. — Orbital phases of photon arrival times. The horizontal lines show the phase 
intervals of the two observations. The shaded areas show the phase intervals with no photons 
detected. The vertical dashed lines show the phases of superior conjunction of pulsar A 
(slightly different in the two observations because of the relativistic advance of periastron). 
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Fig. 3. — Combined light curves for J1537 with three bins. The solid line is the actual counts 
histogram, whereas the dashed lines show the 68% confidence interval for each bin, assuming 
Poissonian statistics. Two periods are shown for clarity. 
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Fig. 4. — Locations in the orbit where the detected photons were emitted: the tickmarks 
on the B orbital circle show the phases of the detected photons. The orbital configuration 
is shown for the two observation epochs, 2009 (top) and 2005 (bottom; as in KPG06). The 
solid vertical lines are the lines of nodes, and the dashed straight lines are the apsidal lines, 
whose orientations are different in 2009 and 2005 because of the relativistic precession. The 
locations of the conjunctions are shown with circles (filled for superior, open for inferior). 
The hatched region of the NS B orbit in the bottom panel shows the section of the orbit 
with no valid data. The shaded sectors show the ranges of true anomalies in which the lack 
of detected photons could be explained by NS B being out of the equatorial wind of pulsar 
A, because of a misalignment of the A's spin axis and the direction of the binary angular 
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Assuming the null hypothesis that the photon arrivals are uniformly distributed in time, 
the probability of detecting n photons within a time interval t is given by the Poissonian dis- 
tribution Pp i SS (n, At), where A is the mean count rate (see Equation (A3) in the Appendix). 
In our case, it can be estimated as A = (Aq + N 2 )/{Ti + T 2 ) = 16.154 counts per binary 
period P or b, where N± = 17 and N 2 = 16 are the numbers of counts detected in the time 
intervals T\ = 0.9925P or b and T 2 = 1.0503P or b, respectively, and P or b = 36,352 s. 

The probability of finding an interval between two neighboring events larger than (or 
equal to) r at a given event rate A is given by Equations (A4) and (Al). For the first 
observation (in which the largest empty interval is Ate = 0.352P or b), these equations give 
the probability Pi(n > 0.352P orb ) = 3.81%. 

For the second observation, the largest intervals, Ati = 0.140P or b and At 17 = 0.134P orb , 
are between the beginning of observation and the first count and between the last count and 
the end of observation. In this case, the exposure time in the gap is the sum of the two 
intervals at the start and end of the observations (see Figure [2]), Ati + Atn = 0.274P or b. 
The probability of obtaining such an interval by chance for uniformly distributed counts is 
estimated in the Appendix: P 2 (r 2 > 0.274P orb ) = 19.0%. 

We must assess the probability of finding this set of two gaps together, i.e., how likely 
is it that we will find by chance gaps as large or larger than observed in two subsequent 
observations if the photons arrivals are uniformly distributed in time? We have calculated 
the probabilities of measuring a gap of a given (or larger) size within a single observation 
(Pi and P 2 ) above. An intuitive cho ice of statistics for calculating the combined probability 
is then Z = P\P 2 (e.g., Fisher 1970 ). (Note that PiP 2 is not itself the probabilityEI) . Since 



Pi and P 2 are both unifo rmly distribute d, the cumulative distribution function is given by 



P{Z < z) = z — zliaz (jSpringerl Il979l . page 102). In our case, z = 0.00 72 and P(Z < 



0.007) = 0.043. Exactly the same probability follows from the recipe by iFisherl (119701 ): 
P(Z < z) = P(x 2 < Xo> v = where the r.h.s. is the probability of observing a value of 
chi-squared larger than xl f° r v degrees of freedom, Xo = ~ 2hiz = 21nFI^ 1 p, and is 
the number of observations (N = 2 in our case). 

Finally, we account for the fact that the largest interval in the first observation and 
the largest interval in the second observation happen to overlap in orbital phase space by a 
significant margin. For two randomly chosen phase intervals A(f)' and A0", the probability 



2 Note that the probability Pi was incorrectly estimated by KPG06 as Pi = P Po i ss (0, NiAte/Ti) = 0.34%. 

3 For instance, for n observations the corresponding product Y[i Pi w iU always tend to zero for large n 
(even if the measured random process samples the null hypothesis) since each probability Pj < 1. 
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of overlap by at least an amount <j) is 

P = A<j>' + A0" - 2<j) , (1) 

for < O < min(A0', A0"). In our case, Ac// = 0.352, A0" = 0.223, and O = 0.177 (see 
Figure 2), leading to the overlap probability P a = 22.1%. 

The combined probability of the two gaps being at least as large as observed and over- 
lapping, is 

P = PcombPo = 0.22 x 0.043 = 0.0094 (2) 

Thus, with a 99.1% probability (2.6<r significance), the observed two large, overlapping gaps 
suggest a physical dependence of the count rate on the orbital phase rather than a random 
fluctuation in the count rate of a steady source. We confirmed the validity of the above 
calculations by extensive Monte-Carlo simulations. 

Finally, we used an alternative method based on the binned light curve (shown in Figure 
|3]) to estimate the probability of the observed count deficit. We produced the phase- folded 
light curve by binning source counts from both observations. Note that we chose the bin with 
the deficit to be centered precisely on the physically significant phase of apastron, (pA = 0.5, 
rather than picking a special phase which minimizes the number of counts in the bin. The 
Poissonian probability of a bin containing only 3 or less photons when expecting 11 is 0.49%. 
We can also evaluate the likelihood of such a light curve being a realization of a constant-rate 
model, with the expectation number of counts in each bin = 33/iV for N phase bins. Since 
the numbers of counts per bin are small, the commonly used \ 2 st atistic (app ropriate for 



large numbers of counts) is biased. Instead we can use the modified ICashl (119791 ) statistic 



£7 = 2 j: 



N 



rii x 
i 



rii \ n. 



In 1 + 



(3) 



where rii is the measured number of counts in the zth bin. For 3 bins, we find C = 10.95. 
To test how likely it is to find a higher value of C by chance, we generate 10000 sets of 
photon arrival phases, drawn from a uniform distribution in the 0-1 phase range. For each 
realization, we generate binned light curve as above, and calculate C. We find that C > 10.95 
in 0.4% of cases. Thus, the null hypothesis of a count rate independent of phase is ruled out 
with 99.6% confidence. The apparently higher significance reflects our choice of the number 
of bins in the light curve. 

To visualize the placement of the observed deficit in relation to the NS orbits and the 
observer's line-of-sight, we plotted in Figure 4 the configurations of the orbits and marked 
the true anomaly values corresponding to the photon arrival phases. The relativistic advance 
of periastron by 5u = uT = 1.756° yr^ 1 x 4.465 yr= 7.84° is accounted for. 



2.2. Spectrum 



We extracted the spectra and produced response files following standard ACIS data 
reduction procedure. Spectral analysis was performed with Sierp^l in the 0.3-8 keV range. 
Due to the low numbers of photons, we fitted the models by minimizing the C-statistic (Cash 
1979) without binning. We fitted a power-law (PL) model, modified by interstellar absorption 
(equivalent to wabs*powerlaw in XSPEC), and fixed the absorbing column to Nh = 3.6 x 
10 20 cm -2 . This value corresponds to the pulsar's dispersion measure for 10% ionization, 
and it is close to the Hi column density through the Galaxy in this direction. Note that 
such a low extinction makes virtually no effect on the spectrum above 0.3 keV. We find that, 
between the two data sets, the best-fit parameters are consistent within their la uncertainties 
(see Tabled]). The combined fit results are shown in Table [U and the confidence contours 
corresponding to these parameters are plotted in Figure O The unabsorbed luminosity is 
£o.3-8keV = (4.0 ± 1.0) x 10 29 ergs _1 , for the distance of 1.02 kpc. 

We also jointly fitted an absorbed blackbody (BB) model and obtained a slightly worse 
fit (C = 154 vs. C = 127 for the PL). The best fit BB parameters, kT = 262^° e V and 
R = 20^g 5 m, are at the hotter /smaller end of the confidence regions calculated by KPG06. 
We cannot exclude the possibility that the spectrum is described by a multi-component 
model, but the statistics are too poor to analyze such models. 



3. Discussion and Conclusions 



KPG06 suggested three possible scenarios to explain the X-ray emission from J1537: (a) 
the emission is generated in the magnetosphere or heated polar caps of the recycled pulsar 
A; (b) the emission results from interaction of the pulsar A's wind with the NS companion, 
such that the strength of the interaction depends on the varying distance between the two 
NSs in the eccentric orbits; and (c) emission only happens when the companion NS B passes 
through the equatorial outflow of the A's wind. Model (a) predicts no dependence on orbital 
phase, (b) predicts fainter emission at apastron and brighter emission at periastron, and (c) 
predicts two segments with faint/no emission at opposite phase intervals, around apastron 
and periastron. Note that modulation would be possible for scenario a) via absorption in the 
magnetosphere of NS B; however with the small magnetosphere of NS B and the inclination 
of the orbit, it is impossible for this system. We find a deficit of photons around apastron, and 
no evidence for enhanced emission at periastron. Thus, our data appear to support model 
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Fig. 5. — Confidence contours (68% and 90%) for absorbed PL model fit to both data 
sets jointly, for photon index T and normalisation Af. Dashed lines show loci of constant 
unabsorbed flux, labelled by flux values in units of 10~ 15 ergs _1 cm~ 2 , for 0.3-8 keV. 



Table 1. Spectral fits 



Parameter 01d a New Joint 



T 2.9±0.5 3.3±0.6 3.07^ 

J\f b 6.4iJ;| 5.7±i:f 6.11};? 

N H (10 20 cm- 2 ) 3.6 3.6 3.6 

^o.a-skev c 3.2±0.9 3.3±1.1 3.2±0.8 

C re d d 0.14 0.10 0.12 



Note. - All fits are performed in the 0.3- 
8keV energy range, the background ignored, and 
C-statistics are used without binning. Nh is held 
fixed. 

a Fit slightly different from KPG06, due to the 
newer callibration used and an extra photon in- 
cluded. 

b Normalization in units of 

10~ 7 photons keV _1 cm _2 s _1 at IkeV. 

c Flux in units of 10~ 15 ergs _1 cm -2 , in 0.3- 
8keV range. 

d C-statistic JCashl Il979h per detector energy 
channel. 
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(c). As the pulsa r's equatorial plane and the orbital plane are misaligned (by 5 = 25° ± 4°; 
Stairs et al.ll2004l ). the NS B will intersect the wind for only part of its orbit, consisting of two 
segments (see Figure 4). One expects more emission as the NS B plunges through the pulsar 
wind, and less when it is outside the wind. Following the calculation in KPG06, the current 
results support a limit on the equatorial outflow opening angle of about A > 20°. X-rays can 
be generated in the intra-binary shock or can originate from the J1537B's magnetosphere or 
surface, bombarded by the particles of J1537A's wind. 

We should note that the phases of superior conjunction of pulsar A (slightly different 
in the two observations because of the relativistic advance of periastron) are also within the 
A(j) w 0.18 phase interval in which no photons were detected in the two observations (see 
Figures 2 and 4). However, since an excess of flux around superior conjunction would be 
more natural than a deficit, due to the Doppler boost, we believe that the observed deficit 
is connected with NS B being outside of the pulsar A's wind around apastron, and is not 
related to superior conjunction. 

In the light of this interpretation, it is illustrative to compare the X-ray properties of 
J 1537 and J0737, the only other DNSB detected in X-rays. For J0737, the pulsed flux 
from the recycled pulsar J0737A accounts for « 76% of th e total DNSB flux in the energy 
range 0.15-4keV, the pulsed fraction growing with energy (jPellizzoni et al.ll2008l ). There is 
also some evidence for pulsations from the slower, non-recycled pulsar J0737B, but only for 
small (A(j) = 0.25) part of the orbit, and the pulsed fraction is rather uncertain. Measuring 
pulsations with the J1537A's period in the J1537 X-ray emission would provide a lower limit 
on the J1537A's contribution to the total observed X-ray flux, but the low time resolution 
of the Chandra ACIS observations does not allow such a measurement. 

In contrast with J1537, no orbital phase depende nce is seen in the X-r ay emission of 
J0737, down to the limit of < 15% on pulsed fraction (jPellizzoni et al.ll2008l ). It itself does 
not exclude a sizeable fraction of the JO 73 7 luminosity to be pow ered by the J0737A's wind. 
Indeed, both the eccentricity (e = 0.0877; iPossenti et al.ll2008l ) and t he angle between the 
orbital angular momentum and the A's spin (< 14° at 95% confidence; iFerdman et al.ll2008l ) 
are much smaller than those in J1537. As a result, the distance between the pulsars changes 
very little, and pulsar J0737B is always near the equatorial plane of J0737A, i.e., within the 
wind of J0737A, even if the opening angle of the equatorial outflow is small. One would not, 
therefore, expect the emission caused by the interaction of the J0737A's wind with J0737B to 
vary significantly with orbital phase (except, possibly, for variations caused by the changing 
view of the shock front). 

The spec tra of both J1537 and JO 73 7 are rather soft in comparison with a typical pulsar 
wind nebula (IKargaltsev fc Pavlovl 120081 ). and their slopes, T « 3.1 and T « 3.4 (for the 
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phase-averaged spectruu@; Pellizzoni et al. 2008 ). respectively, are at the upper end of the 
r ~ 2-3 range typical for recycled and non-recycled old pulsars whose spec tra can often be 
approximated by an abs orbed PL if the number of counts is smal@ (e.g., Bogdanov et al. 
20101 ; iPavlov et all 120091 . and references therein). The similarity of the J1537 and J0737 



spectra seems contradicting the interpretation of the J1537 X-ray emission as powered by the 
J1537A's wind because the X-ray luminosities of the two DNSBs are dominated by physically 
different mecha nisms. However, even thoug h most spectra of P WNe around isolated pulsars 
fit harder PLs (IKargaltsev fc Pavlovl 120081 ) . we know very little about emission from pulsar 
winds shocked in the vicinity of their parent pulsar due to the presence of a binary companion. 
It is possible that spectra of intrabinary PWNe are softer than those of ordinary PWN (e.g., 
due to a smaller amount of the magnetic field energy being transfer red to the particles, 
which also implies a higher magnetization, a; cf. [Pellizzoni et al.1 120081 ). We also note that 
the J1537's spectrum is marginally harder than that of J0737, and the large uncertainties 
leave room for even more appreciable differences. Deeper X-ray observations of J1537 should 
better constrain its spectrum and orbital dependence, 

A significant difference between the two DNSBs is that the X-ray radiative efficiency, 
%3-8keV = -^0.3-8 kev/-^A ~ 2.2 x 10~ 4 for J1537, is a facto r of 7 lower than for J0737 
(with account for the revised distance to J0737, d ~ 1150pc; iDeller et al.1 120091 ) . The X- 
ray efficiency of J0737, ?7n. 3-8kev ~ 1-6 x 10~ 3 , is typical for millisecond pulsars detected 



in X-rays (see Figure 4 in IPavlov et al.1 120071 ). in accordance with the conclusion that the 
J0737's emission is dominated by pulsar J0737A (IPellizzoni et al.l 120081 ). On the contrary, 
the J1537's efficiency is lower than those of nearly all other recycled pulsars. If most of the 
X-ray luminosity in J1537 comes from the interaction of the pulsar A's wind with pulsar B, 
as follows from our interpretation, the efficiency of J1537A is even lower, making it one of 
the least efficient recycled pulsars known. 

An upper limit on the wind-powered J0737's luminosity, ~ 1.7 x 10 30 erg s -1 in the 
0.3-8 keV band (which includes the J0737B's luminosity), can be estimated as the difference 
between the total observed luminosity and the luminosity pulsed with the A's period. This 
limit is a factor of 4 larger than the total luminosity of J1537, but this fact does not exclude 
the possibility that the unpulsed luminosity of J0737 is wind-powered because the J0737 



IPellizzoni et al.1 (|2008l ) also performed phase-resolved spectroscopy of J0737A and concluded that three 
components, two BBs (with kT ~ 150 eV and kT = 30 — 85 eV) and a PL (with Y — 3.3^$ 2)) are required 
to simultaneously fit the spectra in three phase intervals. The modulation of the PL normalization was by 
itself able to explain the large pulsed fraction (75.7% ± 5.4% in the 0.15-4 keV band) at J0737A's period. 

6 The spectra of ordinary old and recycled pulsars with high S/N usually require three components for an 
acceptable fit, two thermal components (BB or NS atmosphere) and a PL. 
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orbit is tighter. 

In summary, we have confirmed the deficit of X-ray emission in the J1537 DNSB around 
apastron. This implies that a significant fraction of the X-ray emission is powered by the 
equatorial wind of the recycled pulsar J1537A, and the emission probably occurs when the 
companion NS crosses this wind. Our interpretation also implies that the apparent X-ray 
efficiency of pulsar J1537A is very low compared to other millisecond pulsars. 

We are deeply indebted to Ingrid Stairs for providing us with the up-to-date radio 
ephemeris. This work was supported by Chandra award GO9-0073X and NASA grant 
NNX09AC84G. The research was also partly supported by the National Science Founda- 
tion under Grants No. 0908733 and 0908611. The work by GGP was partly supported by 
the Ministry of Education and Science of Russian Federation (Contract No. 11.G34.31.0001). 
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A. Probability calculation 

The probability that, for an observation with N counts detected, the largest of the N + l 
time intervals between neighboring events (which include the beginning and the end of the 
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observation) is as large as or larger than r can be calculated as the difference between unity 
and the probability of all the intervals being smaller than r: 

P(At >t\N;T) = 

1 - N\T~ N jf d(Atx) ... d{At N ) ®(t-J2 e^AU + T-T^J, (Al) 

where T is the duration of the observation, 6(x) is the Heaviside step-function (0(x > 
0) = 1, Q(x < 0) = 0), and At; is the interval between the (i — l)th and iih events (i — 
and % = N + 1 correspond to the beginning and the end of the observation, respectively). 
Equation (Al) takes into account that, for N fixed, each possible configuration of intervals 
ought to have the same probability of occurrence. The step-functions in Equation (Al) 
correspond to the following constraint 



iV 



< T- ^AU = At N+l < t. (A2) 
i=i 

The analytical calculation of the A^-dimensional integral in (Al) leads to very cumbersome 
expressions; in particular, for each given N, they are different for the separate intervals 
(k + 1) _1 < t/T < fc -1 , where k = 1,2,...N. However, the integral can be evaluated 
numerically. 

The number of counts is itself a random quantity that obeys some statistical distribution. 
For the null hypothesis we are investigating, the probability of detecting n photons in a time 
interval t is given by the Poissonian distribution 

(Xt) n e~ xt 

P Poiss (n,Xt) = , (A3) 

71/. 

where A is the mean count rate. For a given A, the probability that the largest interval 
between neighboring events is as large as or larger than r is the expectation value of P(At > 
t\N;T): 

oo 

P(At >r\\T) = J2 P ( At >r\N;T) P Poiss (iV, AT) (A4) 

N=0 

For the first observation (r = 0.352, T = 0.9925, A = 16.15; the unit of time is the orbital 
period) this gives P(At > 0.352 | 16.03) = Pi = 0.0381. 

For the second observation, the largest interval in the phase space is not the one between 
the consecutive events but the sum of the waiting times before the first count and after the 
last one (see Figure 2). The probability of the sum of two intervals to be larger than r when 
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the total number of events N in the observation time T is given can be calculated as 

P(h + t N+1 >r\N;T) = 
1 - N\T- N d(Ah) d(At N ) (t - ^ ( T " X> - r ) > ( A5 ) 

Since the expectation value of P(ti+tN+i > r | N; T) is P(ti+t N+ i > r \ XT) = J2n=o P(ti + 
*;v+i > r | iV; T)P Poiss (iV, AT), for the second observation we obtain P(ti+t 17 > 0.274 | 16.70) = 
p 2 = 19.0%. 



